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Abstract 

Under the assumption of a smooth full phase-space distribution function we prove that the 
nonequilibrium entropy S which is considered as a functional of the distribution vector for an 
N-body system possesses a lower bound and therefore can not decrease. We also compute the 
rate of change of S , dS/dt, showing that this is non-negative and having a global minimum at 
equilibrium. As an aplication we obtain a generalization of the (Bhatnager-Gross-Krook) BGK 
relaxation model. 
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I. INTRODUCTION 



As it is known, the second law of Thermodynamics has been based on Kelvin's and 
Claussius's principles which establish that work can not be produced without spending 
energy. These principles are based on experiments but to my knowledge no theoretical 
derivation based on first principles has been given^J. For an adiabatically isolated system the 
second law establishes that the entropy can not decrease and is maximum at equilibrium 3]. 

To provide a theoretical proof of the second law or the law of increase of entropy for 
an isolated system first it is necessary to properly define the nonequilibrium entropy which 
would coincide with the thermodynamic entropy at equilibrium. 

The equilibrium entropy is a well established concept. According to Boltzmann the 
entropy of a microstate X of a macroscopic system is a number, Sb(X) = /c^log |Fm(x)|, 
with \T M \ being the volume of the region of the phase space T M corresponding to the 
macrostate M = M(X) and fee is the Boltzmann constant. The macrostate M is all of 
a group of states Y such that M(Y) = M(X) = M. For Boltzmann, every adiabatically 
isolated system in its natural evolution moves towards an state of maximum disorder -the 
state of equilibrium, which corresponds with the greatest possible volume in phase space. 
Gibbs, on the other hand, defines entropy as a functional of the distribution function F 
in phase space, Sq = —k B J FliaFdX, and both entropies coincide at equilibrium. Out of 
equilibrium the situation is more ambiguous, due to the fact that the Gibbs entropy is a 
constant of motion under the microscopic dynamic given through the Liouville equation, it 
has been thought that macroscopic processes require coarse-graining. This point was worked 
out first by Gibbs and P. Ehrenfest & T. Ehrenfest. On the other hand, some authors argue 
that the relevant entropy for understanding thermodynamic irreversibility is the Boltzmann 
entropy and not the Gibbs entropy 

In this scenario, unlike treatments based on the definition of coarse-grained entropiesj 
here we propose a generalization of the fine-grained Gibbs entropy postulate in terms of the 
set of reduced distribution functions of the system, the distribution vector, as the nonequi- 
librium entropy of an N-body isolated system Q], Under the assumption of smooth 
phase-space densities we prove that this nonequilibrium entropy can not decrease and is not 
a constant of motion under the dynamic given through the BBGKY (Bogolyubov-Born- 
Green-Kirkwood-Yvon) hierarchy of equations [9]. Here, we assume smoothness to exclude 
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multifractal phase-space distributions which some simulations have shown|7|. 

It is a fact that the Liouville equation describes an incompressible flow of representative 
points in phase space, however this is not the case with respect to the BBGKY hierarchy. 
Additionally, unlike an incompressible flow a compressible flow produces entropy, thus in this 
fact resides the reason of the behavior of the entropy we propose . Moreover, this entropy 
coincides with the Thermodynamic entropy at equilibrium. Hence, we are able to compute 
the entropy production for an isolated system, which is non- negative (zero at equilibrium), 
therefore conferring the character of the nonequilibrium potential of the system to this 
entropy production. Is this a theoretical proof of the second law of Thermodynamics? it 
depends to the credit one gives to the nonequilibrium entropy we propose, nonetheless this 
is the closest definition to the microscopic dynamics of the isolated system we know. 

We have structured the paper as follows. In section 2, we introduce the Hamiltonian 
dynamics of the N-body system obtaining the generalized Liouville equation. In section 3, 
we analyze the properties of the nonequilibrium entropy, compute the entropy production 
and obtain the kinetic equation for the one-particle reduced distribution function. Finally 
in section 4, we emphasize our main conclusions. 



II. HAMILTONIAN DYNAMICS 



We will consider a N-body system whose Hamiltonian results from the addition of a 
kinetic energy term and a term coming from the interaction between the particles 

N 2 1 N 
3=1 &k=l 

with m being the mass of a particle, and <f>(\qj — q^l) = <pjk the interaction potential. 
Moreover, the equations of motion are 

dH dH 

The statistical description of the system can be performed in terms of the full phase-space 

distribution function F(x N ,t), where x N = {x\, ,xn} and Xj = (qj,Pj) or alternatively 

in terms of the distribution vector f [9], both completely equivalent but the last more 
appropriated for nonequilibrium systems in order to show irreversibility. Here, 

f={fo fi(x u t), f 2 (x 2 , t), f N (x N , t)} (3) 
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is the set of all the s-particle reduced distribution functions, s = 0, ,N, where the s- 

particle reduced distribution functions 

$ s = (jv~7ji / F ( xN ^) dx s+1 ...dx N , (4) 

are obtained by integrating over N — s particles with f — 1. The dynamics of the reduced 
distribution functions can by obtained from the Liouville equation 
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-F=[H,F} p , (5) 



where f..,..L is the Poisson bracket, by using Eq. JU). One obtains the generalized Liouville 

D □ 

equation [7|]- [9( 

|f(*)=£f(t) , (6) 

which is a compact way of writing the BBGKY hierarchy of equations. In this equation, £ 
is the generalized Liouvillian whose diagonal part VC is defined through j^] 

(s\v£\ s >) = 5 s ,Jj2 L ° + E M > ( ? ) 

I j=l j<n=l ) 

where \s) represent the s-particle state defined through (s| f = f s (x s ), with x s = 

{xi, ,x s }, s = 0, ,N. The s-particle states have the property (s \ s') = 5 StS i such 

that J2^=o I s )( s \ = 1 constitutes the closure relation. Moreover, L° = [fZJ, ...] , where 
H° = p'j/2m, and L'j n = [H'^ n , ...] p , with Hj n = \<\>^ n - In addition, the non-diagonal part 
QC is given by 



(s\Q£\s') =5 s , tS+1 j |E L j>+i| dx ^ ■ 



(8) 



Finally, in view of Eqs. (|IJ) and (jHJ) and using the properties of the s-particle states \s), by 
multiplying by \x s ) to the right and by (s\ to the left of © one obtains one of the components 
of the hierarchy 

£~\ S S S n 

-f s -J2 L °jf° - E L 'U< E / L'^+ifs+idx^ , (9) 

j=l j<n=l j=l 

where the right hand side is the collision term coupling f s with / s +i. In the case of s = 1 
Eq. (jUJ) reduces to 



where F x 2 = — V 



1*1,20- 



p i+Pl l- h = -J Fh2 A- f2dX2 , (10) 



III. ENTROPY PRODUCTION AND IRREVERSIBILITY 



Our statistical theory is based on the definition of the nonequilibrium entropy proposed 
in fl, 3 



S = -k B tr {fin (f^f)}+S eq 

= -^b^— / /shiy — dxi dx s + S eq (11) 

s=1 s - J Jeq,s 

which generalizes the Gibbs entropy postulate, where S eq is the equilibrium entropy coincid- 
ing with the thermodynamic entropy and f e(? is the equilibrium distribution vector satisfying 
£f eq = 0. The nonequilibrium entropy (fTTJ) reaches its maximum value S = S eq at the 
equilibrium state, when f = f eq . 

At this point, it should be clarified that a description of the system in terms of the full 
phase-space N-particle distribution function is justified for equilibrium systems for which 
the representative ensemble is distributed uniformly through the phase space. Nonetheless, 
the nonuniformity of nonequilibrium systems leads to a random clusterization caused by the 
interaction between different s-particle clusters, and the same global quantity may be ob- 
tained for an infinite number of realizations. This fact is taken into account in the BBGKY 
hierarchy jllT] which explicitly incorporates the interaction between different s-particle clus- 
ters through the nondiagonal part QC of the generalized Liouvillian C The nonequilibrium 
system behaves as a random mixture of s-particle systems and, therefore the total entropy 
should be the sum of all the s-particle entropies. The randomness creates entropy as we will 
see later on. 

More interestingly here is to establish the direction of change of this entropy in a natural 
process. To discern this we first note that since f s In f s /f eq , s is a convex function[12] 

f a In £ f s - f eqjS ( for f s ^ and f eq , s > 0) (12) 

Jeq,s 

and due to the fact that f s and f eq<s both are smooth phase-space densities in the s-particle 
phase space 

// s ln— — dx\ dx s ^0. (13) 
Jeq,s 

Thus, given that f s is bounded and since f eq ^ s > there must exist a constant C such that 

/ s hiy — dx 1 dx s < C , (14) 

Jeq,s 



for s = 1, N. Hence, a lower bound for the nonequilibrium entropy exists 

N 1 

0>S-S eq ^-k B J2~C , (15) 
s =i s - 

which means that the nonequilibrium entropy (fTTj) can not decrease. Consequently, the law 
of increase of entropy (i.e. the second law of thermodynamics) follows, constituting this 
result one of the goals of the present contribution. 

The rate of change of the nonequilibrium entropy (|11|) or entropy production can be 
computed using Eqs. (JHJ), giving 
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= -^r{|ln(f-/f)} = 

-Afetr^flnft^f)} , (16) 



which according to Eq. (JT5|l should be non-negative. In a more explicit way, after using Eq. 
@, Eq- (HH can be written as 

d S 1 X ^ f ( d s i 

-qt = -f Yl 7\ Yl / /*pj ~ kBT 7fo. ln + 2 F ^ + ^ dXs 

s=l ' j=i ^ \ Hj' jVi=l / 

where J-j(x s ) is defined through f s (x s )J z j(x s ) = J Fj . s +if s +\dx s+ i, T is the equilibrium 
kinetic temperature taking into account that the dependence of / e?)S in the velocities is 
given through a local Maxwellian and F^j = —^j4>j,i- 

The entropy production given through Eq. (|T7|) vanishes at equilibrium when f s coincides 
with f eqyS . Moreover, because pj is arbitrary 

is sufficient to satisfy the extremum condition 5SJ6f s \ eq — 0, with 5 = dS/dt. Here, the 
right hand side of Eq. (fT%|) is the mean force [l2j and the left hand side is the sum of 
two parts, the first is the force due to the van der Vaals interactions with the s — 1 fixed 
particles different from the j-th particle in the s-th cluster (responsible for the compressions 
and dilatations of the s-th cluster), while the second is the average force on particle j-th 
by the remaining N — s particles of the system. It is known that Eq. (JTSj) gives rise to a 
hierarchy of equations the Yvon-Born-Green (YBG) hierarchy [21]. Now, by making use of 
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the extremum condition Eq.(JTHJ), we can rewrite the entropy production Eq. (jl7j) as 



^ = E / /'Pi fa' ~ T ?) dxi dXs 



3=1 
N s 

- f E 7? E / ( A /« + /«*-) p^ fa - ^T) > ( 19 ) 

s=l ' j=l J 

where A/ s — f s — f eq , s - The right hand side of the last equal sign in Eq. (fT9~)) has two 
parts, the second of which involving f eq>n , provided T 3 - has even parity under time reversal, 
vanishes. It must be pointed out that Tj has even parity if f s+ i(x s+1 ) = fi(x s+ i)f s (x s ) 
which is a reasonable assumption in large enough time scales. Henceforth, after defining 
the nonequilibrium current Jj = A/ S pj which vanishes at equilibrium, Eq. (JT9*)l can be 
rewritten 

§ - E ji E / *i fa - *T) • ( 20 ) 

8=1 j=l 

In the previous equation (|2*0*j) the term under the integral possesses the form of the prod- 
uct of a thermodynamic current times the corresponding conjugated thermodynamic 
force [T d - — ^j 9 ) , which constitutes the typical expression of the entropy production in 
Nonequilibrium Ther mody namics. At this point, following the standards of Nonequilib- 



rium Thermodynamics 
assumed, it follows 



14j . from Eq. (|20|) provided the existence of an equilibrium state is 



-E^fa-^f), ( 21 ) 



T 

3=1 



where is a phenomenological coefficient which should be non-negative and having even 
parity or in terms of the related mobility Mjj = ~Lij/Tf s 

s 

Ji = - E ^ fa - ^ ■ ( 22 ) 
i=i 

The phenomenological relations 1)21)1 or ()2~2~|) manifest the fact that a restoring current arises 
when the system deviates from equilibrium by the action of a thermodynamic force, this 
constituting a generalization of the Onsager's regression hvpotesisjl^]. Additionally, My 
might be a function of f s and {Fj}, Mjj(/ S , {^}}) in general. Therefore, after substituting 
Eq. (|22D into Eq. (JHD 

§ = f E ^ E / fa - ) 2 ^ • (23) 

8=1 ' i,i=l ^ 



To test our previous analysis, near equilibrium we can define an effective mobility = 
/ Mijfeg^dxi dx s , thus 

j, -^./;m;;(^-^;") (24) 

3=1 

which by introducing the inverse mobility matrix (Yl^iCij^-jf = ^u)i the friction ma- 
trix, can be rewritten as 

3=1 

-Kft^-W • (25) 

3=1 

For s — 1, the previous equation allows us to derive from Eq. (fTU|) the following kinetic 
equation for f\ = f 

y + p ^ + kBT (*. ln /ea ) b_, = c ^p (/ _ M (26) 

constituting a generalization of the BGK [15j relaxation model containing the mean force 
and an extra acceleration term involving d/dp(f — f eq ) which distinguishes this equation 
with respect to the original BGK. 

IV. CONCLUSIONS 

In the framework of the BBGKY description we have found a way of generalizing the 
Gibbs fine-grained entropy for an isolated N-body system leading to irreversibility. We have 
proved that the nonequilibrium entropy of an isolated N-body system defined as a convex 
functional of the distribution vector can not decrease and that the rate of change of this 
entropy is a non-negative quantity having a global minimum at equilibrium. As an appli- 
cation of those results, after computing the entropy production, following the methods of 
Nonequilibrium Thermodynamics we introduce a positively defined mobility matrix which 
relates the nonequilibrium current with the corresponding thermodynamic force enabling us 
to write the entropy production as a positively defined quadratic form. As a consequence 
of introducing the phenomenological law, we have obtained a generalization of the BGK 
relaxation model which contains the mean force and an extra acceleration term involving 
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9/ dp (/ — f eq ). At this point it should be mentioned that the nonequilibrium thermodynam- 
ics analysis may be generalized to the phase space as done in the Mesoscopic Nonequilibrium 
Thermodynamics (MNET) 17 1. 

It should be stressed that as we have shown a proper description of the nonequilibrium 
system might be given in terms of the distribution vector which evolves according to the 
genreralized Liouville equation. The compressibility inherent to the BBGKY hierarchy leads 
to the randomness, so that the system can be thought either as a random mixture of s- 
particle systems or an interacting mixture of compressible fluids. Hence, this description 
brings about the correct macroscopic dynamics. 

Due to the generality of our theory of irreversibility even the Boltzmann irreversible equa- 
tion and hence the Boltzmann H-theorem are embedded in here, as was shown in a previous 
paper 0. Moreover, a quantum formulation is possible, which will be shown elsewhere. 
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